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Abstract 

The study of the formation of molecular hydrogen on low temperature surfaces is of interest 
both because it enables to explore elementary steps in the heterogeneous catalysis of a simple 
molecule and due to its applications in astrochemistry. Here we report results of experiments 
of molecular hydrogen formation on amorphous silicate surfaces using temperature-programmed 
desorption (TPD). In these experiments, beams of H and D atoms are irradiated on the surface 
of an amorphous silicate sample. The desorption rate of HD molecules is monitored using a mass 
spectrometer during a subsequent TPD run. The results are analyzed using rate equations and 
the energy barriers of the processes leading to molecular hydrogen formation are obtained from 
the TPD data. We show that a model based on a single isotope provides the correct results for 
the activation energies for diffusion and desorption of H atoms. These results are used in order to 
evaluate the formation rate of H2 on dust grains under the actual conditions present in interstellar 
clouds. It is found that under typical conditions in diffuse interstellar clouds, amorphous silicate 
grains are efficient catalysts of H2 formation when the grain temperatures are between 9-14K. This 
temperature window is within the typical range of grain temperatures in diffuse clouds. It is thus 
concluded that amorphous silicates are good candidates to be efficient catalysts of H2 formation 
in diffuse clouds. 
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I. INTRODUCTION 



Few are the studies of the formation of molecular hydrogen on low temperature sur- 
faces. One of the pioneering experiments was done in the 1970's by the group of Giacinto 
Scoles, who measured the scattering, sticking and energy deposition of atomic and molecular 
hydrogen beams on the surface of bolometers (semiconductor thin films) at liquid helium 
temperature It was found that both the sticking coefficient and the hydrogen 

recombination rate depend on the coverage of H2 on the target surface. It was also shown 
that the heat released in the formation of molecular hydrogen causes the desorption of hy- 
drogen molecules that have been pre-adsorbed on the surface. Thus, a molecule just formed 
is immediately ejected from the surface. These experiments offered a rare view of the inter- 
action of hydrogen atoms and molecules in the physical adsorption regime and a connection 
with processes in interstellar space. However, the sample temperature of 3-4 (K) was well 
below that of interstellar dust grains, the coverage of the sample with atoms/molecules was 
high and the ice layer not fully characterized. These conditions made it difficult to obtain a 
quantitative understanding of actual gas-dust grain processes in astrophysical environments. 

Molecular hydrogen (H 2 ), the most abundant molecule in the Universe, influences the 
chemical make-up of the Cosmos 4, [5| and is instrumental in the formation of stars by con- 
tributing to the cooling during the gravitational collapse of molecular clouds. The challenge 
of explaining the formation of molecular hydrogen in space begins with the realization that 
the stabilization of the nascent molecule in the bonding of two (neutral) hydrogen (H) atoms 
involves the forbidden transition to the ground state. Three-body gas-phase interactions are 
too rare to contribute significantly to H 2 formation in cold clouds 4j, but may take place 
in other environments such as interstellar shocks. Under conditions observed in interstellar 
clouds, other gas-phase routes (such as: H+e — > H + hu, H + H — > H 2 + e; or, less 
frequently, H + H+ -. H+, H+ + H ^ H 2 + H+) do not tnake enough H 2 to counterbalance 
the known destruction rate due to UV photons. |4J. 

In the 1960's, Salpeter and collaborators proposed a model in which H 2 formation occurs 



on the surfaces of interstellar dust grains 
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,|9j. These grains are formed in the envelopes 



of massive late-stars and in novae and supernova explosions. They are made of carbonaceous 
materials and of silicates. Their sizes exhibit a broad power- law like distribution between 
1-100 nm [lO, 



11 



121 ] . Observations of scattered, absorbed and emitted starlight, and lab- 
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oratory work, show that in the interstellar medium (ISM), silicate grains are amorphous 
and mostly of composition (Fe :r Mgi_ x ) 2 Si04, where < x < 1 [131 ] . There is, on average, 
one dust grain per about 10 12 hydrogen atoms, and the grains account for about 1% of the 
mass of interstellar clouds. Kinematic calculations show that in order to produce enough 
molecular hydrogen to counterbalance the destruction rate, the catalysis on grain surfaces 
must be efficient. More specifically, the processes of hydrogen sticking, migration and bond 
formation on the grains must convert at least about ~ 30% of the adsorbed hydrogen atoms 
into molecular form [3]. 

With some exceptions Q, [l^ , chemical models that look at the chemical evolution of an 
interstellar cloud, have largely ignored or underplayed the coupling of gas and dust. However, 
observations, experiments and calculations are pointing to the fact that the formation of 
key ISM molecules [such as H 2 , formaldehyde (H 2 CO), and methanol (CH 3 OH)] takes place 
on dust grains, as gas-phase reactions are too slow in these particular cases 0, Q, Q, 



191 ]. As far as the formation of molecular formation is concerned, there is a great need to 



know the basic mechanisms of reaction (Langmuir-Hinshelwood, Eley-Rideal or hot atom), 
characteristic energies for various processes (diffusion and desorption) and kinetic parameters 
of dust-catalyzed reactions so they can be used in models of interstellar chemistry. 

It is within this framework that in the late 1990's we began a series of investigations on 



21 



221 ] . These experi- 



the formation of molecular hydrogen on analogues of dust grains [2( 
ments, inspired by Scoles' work, were aimed at combining tools of surface science, chemical 
physics and low temperature physics in order to recreate the environmental conditions of the 
interstellar space and overcome some of the limitations of prior experiments, such as: high 
fluxes of H, too low sample temperatures and not adequately characterized materials (For 



a review of early experiments, see Ref. 



23j). In practice, the experiments have to be done 



at low background pressure, low sample temperatures and low fluxes of atoms impinging on 
the samples. The first two requirements are relatively easily achieved. Even taking special 
care to obtain fluxes of low energy (200-300K) hydrogen atoms, it is not possible to either 
produce or detect as low fluxes of atoms as appear in the ISM. Thus, carefully designed the- 
oretical and computational tools need to be used to simulate the actual processes occurring 
is the ISM using the results of the experiments. 

The formation of molecular hydrogen on surfaces has been explored at length in the past, 
but most of the work has been on characterized surfaces of metals and semiconductors, and 
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at much higher surface temperatures and fluxes (or coverages) than in the regime we are 
interested in. On low temperature surfaces, efficient recombination can occur only if the 
mobility of hydrogen is high. There are situations when this does not have to be verified, as 
in the Eley-Rideal and hot atom mechanisms, in which H atoms from the gas phase directly 
interact with the target hydrogen atoms or move on the surface at superthermal energy. 
Such mechanisms have been shown to be working in the interaction of H with H-plated 



metal 
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25| ] , silicon 



26( and graphite 



271 ] surfaces. Although there are certain interstellar 



environments where these mechanisms enter into play, the diffuse cloud environment - where 
the coverage of H atoms on a grain at any given time is very low - is not one of them. Thus, 
we expect that H atoms in our experiments to experience physical adsorption forces and the 
dominant mechanism of reaction is expected to be the Langmuir-Hinshelwood reaction. 

In the experiments, the sample is exposed to well collimated beams of hydrogen (H) 
and deuterium (D) atoms. The production of HD molecules occurring on the surface of a 
dust grain analogue is measured both during the irradiation with the beams and during a 
subsequent temperature programmed desorption (TPD) experiment. In order to disentangle 
the process of diffusion from the one of desorption, additional experiments are carried out 
in which molecular species are irradiated on the sample and then are induced to desorb. 
We first studied the formation of hydrogen deuteride on a telluric polycrystalline sample 



of olivine 



3, 



2l\ . This was 



amorphous water ice 
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bllowed by studies of HD formation on amorphous carbon and 
2i|. Water ice is known [together with other condensables, such 
as carbon dioxide (CO2), carbon monoxide (CO), and CH3OH] to coat grains in molecular 
clouds, under conditions of high density, low temperature (10-15K) and in the absence of 
UV radiation. The high density is important for effective mantle formation, while the low 
temperature is essential in order to avoid mantle evaporation. Thus, shielding is required 
from both UV radiation (which causes photodissociation) and from thermal heating sources, 
namely stars and protostars. Other groups studied the formation of molecular hydrogen 
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on amorphous water ice surfaces [30] as well as the desorption of H 2 from these surfaces 



Measurements of the kinetic energy of hydrogen molecules emerging from amorphous 



water ice show that the molecules have nearly thermal energies 29|, ]30| . The ro- vibrational 



states of excitation o 



graphite samples [30 



the just-formed molecules leaving the surface were also studied using 



34j. For a review of recent experimental work on the formation of 
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molecules on astrophysically relevant surfaces see Ref. [35]. 

The picture that emerges from the experimental studies on H2 formation is the following. 
Molecular hydrogen formation takes place via the Langmuir-Hinshelwood mechanism, at 



least in the range of surface temperatures (during irradiation) that was sampled 



(5-16 



K). The 



results of our TPD experiments were analyzed using rate equation models |36|, |37j, |38|, 139 ]. 
In this analysis the parameters of the rate equations were fitted to the experimental TPD 
curves. These parameters include the energy barriers for atomic hydrogen diffusion and 
desorption and the energy barrier for molecular hydrogen desorption. Using the values of the 
parameters that fit best the experimental results, the efficiency of hydrogen recombination 
on the polycrystalline olivine, amorphous carbon and water ice surfaces was calculated for 
interstellar conditions. By varying the temperature and flux over the astrophysically relevant 
range, the domain in which there is non- negligible recombination efficiency was identified. 
It was found that the recombination efficiency is highly temperature dependent. For each of 
the samples, there is a narrow window of high efficiency along the temperature axis, which 
slowly shifts to higher temperatures as the flux is increased. Contrary to expectations, for 
astrophysically relevant fluxes, the formation of H2 on polycrystalline olivine occurs in a 
temperature range which is too low and too narrow for making molecular hydrogen but in 
very selected environments. On other surfaces, such as amorphous carbon and amorphous 
water ice, efficient molecular hydrogen formation occurs on a wider sample temperature 
ran ge, which is within the typical range of grain temperatures in diffuse and dense clouds 
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38, 31 



Since silicates make up a significant fraction of the interstellar dust, and given our results 
that pointed at a less than optimal efficiency of polycrystalline olivine in catalyzing the 
formation of molecular hydrogen, we decided to revisit the formation of H2 on such a class 
of materials. Although crystalline silicates are observed in circumstellar envelopes, in the 
ISM silicates are mostly amorphous. We thus studied the formation of molecular hydrogen on 
amorphous silicate samples, (Fe^, Mgi_ :r ) 2 Si04, x = 0.5, produced by laser ablation (with 
wavelength of 266 nm) of a mixed MgO, FeO and Si02 target in an oxygen atmosphere 
(10 mbar). The optical and stochiometric characterization of the samples produced with 



this technique is given elsewhere 



40|. Measurements of the amorphous silicate samples by 



scanning electron microscope show rough surfaces with droplets of different sizes. However, 
unlike the amorphous water ice, no indications of porosity were observed. A brief report on 
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the formation of molecular hydrogen on an amorphous sample of composition (Feo.5Mgo.5)2 
SiC>4 at low irradiation temperature (around 5K) appears in Ref. 41]. 

In this paper we present data on the formation of HD at a higher irradiation temperature 
(10 K) which is at the lower edge of the range of the dust temperature in the ISM. We describe 
in detail the methodology for analyzing the results of the TPD experiments. We show 
that the results, based on measurements of HD formation are applicable for the evaluation 
of H2 formation rate in interstellar clouds. Results from surfaces of samples with other 
compositions (Fe^Mg^^SiO/i, < x < 1, will be given elsewhere. 

The paper is organized as follows: In Sec. [TT] we present a review of the experimental 
methods. TPD results for HD formation on an amorphous silicate surface are given in Sec. 
IIHI The rate equation model is presented in Sec. [IV] and is used in Sec. |V]for the analysis 
of the experimental results. Applications to interstellar chemistry are considered in Sec. IVIl 



II. REVIEW OF EXPERIMENTAL METHODS 



The apparatus consists of two atomic/molecular beam lines aimed at a target located in a 
ultra-high vacuum (UHV) chamber. The triple differentially pumped beam lines have each a 
radio-frequency (RF) dissociation source; typical dissociation rates are in the 75-90% range 
and are measured when the beams enter the sample chamber. The atoms can be cooled by 
passing them through a short aluminum nozzle connected to a liquid nitrogen reservoir via 
copper braids. The fluxes F H and F D , of H and D atoms, respectively, are both estimated 
to be equal to F = 10 12 (atoms cm" 2 s _1 ) [28 ] . 

The detector, located in the main UHV chamber, is a differentially pumped quadrupole 
mass spectrometer that can be rotated around the sample; it is used to measure the signals 
proportional to the number of particles in the beams and of the molecules evolving from the 
surface either during the irradiation phase or during the thermal desorption. The sample is 
mounted on a liquid-helium cooled sample holder with thermal shields. The sample can be 
rotated around its axis. A heater in the back of the sample can heat it to 400K for cleaning. 

A typical experiment proceeds as follows. During the baking of the apparatus and prior 
each series of data taking, the sample temperature is taken to ~ 400K. After cooldown, 
when the desired sample temperature is reached, the sample is exposed to two converging 
beams of H and D (prior the pumping down, laser beams are shone through the beam lines 
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to make sure they point at the same spot on the sample). Using two isotopes is essential in 
order that the fraction of undissociated molecules and the background pressure of H2 will 
not affect our measurements. During the irradiation phase of the experiment, as well as the 
subsequent TPD phase, the detector monitors the increase of HD partial pressure due to 
HD formation on and release from the surface. After irradiation time of to (s), the beams 
are turned off and the sample temperature is raised, either by shutting off the flow of liquid 
helium or by increasing the power to the heater. The sample temperature vs. time, T(t), 
is measured by an iron-gold/chromel thermocouple and a calibrated silicon diode placed in 
contact with the sample. The temperature curves are nonlinear but highly reproducible. 
To account for this nonlinearity, the TPD curves in Figs. [20 show the instantaneous HD 
desorption rate vs. temperature rather than vs. time. The heating rate is steep in the 
beginning and gradually decreases. The temperature curves for the experiments analyzed in 
this paper are shown in the insets in Figs. [1] and [2j The typical irradiation times are 1 to 2 
minutes (although, the actual time the sample is exposed to the beam is half that, since the 
beams are chopped with a mechanical selector of 50% duty cycle). The typical coverages, 
assuming a sticking coefficient of 1 are of the order of a few percents of a layer. 

During the experiment, it is checked that after a thermal desorption there is no re- 
adsorption of HD on the surface. Also, in the experiments done at the lowest sample 
temperatures during irradiation (such as 5.6K), care is taken to make sure that the sample 
remains in this low temperature long enough to reach thermal equilibration (as judged by the 
repeatability and quality of the results). Further details can be found in other publications 
42]. 



III. EXPERIMENTAL RESULTS 



Here we report TPD results for HD desorption from amorphous silicates of composition 
(Fe x ,Mg 1 _ x ) 2 Si04, where x = 0.5, after irradiation by H+D beams at a sample temperature 
of 10K. Three specimen were used, yielding comparable results. These results complement 
and extend those reported in Ref . [41] , where we analyzed TPD results for irradiation of H 
and D (hereafter H+D) at a low sample temperature (~5K) and in which the exposure was 
varied. In these experiments we find that most of the molecules that form on the surface 
come off during the TPD; therefore, we focus our analysis on this aspect of the experiment. 
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In a separate set of experiments, beams of HD molecules were irradiated on the same surface. 

In Fig. [1] we show a typical TPD trace of HD desorption, taken after irradiating the 
sample, while at 5.6K, with HD molecules (x) and with H and D atoms (circles). Both curves 
exhibit a strong peak at about 15-16K and a lower but broad peak (or a shoulder) around 
21K. The broad peaks reflect the fact that the surface is disordered. For comparison, a TPD 
trace taken from a polycrystalline olivine sample is much narrower and has a maximum at 



a lower temperature 



4l|. 



In Fig. [2] we present the desorption rates of HD after irradiation with HD molecules (x) 
and H+D atoms (circles) at a sample temperature of ~ 10K. These TPD traces exhibit two 
peaks. The main peak coincides with the high temperature peak of Fig. [TJ In addition, 
there is a smaller peak or a shoulder at ~31-32K, which does not appear in Fig. [TJ However, 
it turns out that the TPD curve of Fig. 1 was recorded only up to 24K. Therefore, we cannot 
reject the possibility that a third peak, around 31-32K also appears under the conditions of 
Fig. 1. 

The experiments with HD irradiation provide direct information on the energy barriers 
for the desorption of HD molecules. Combining the two sets of experiments, we find that 
there are three types of adsorption sites for molecules: shallow, medium and deep. When 
HD molecules are irradiated at 5.6K they reside in the shallow and medium depth sites; 
when the molecules are irradiated at 10K they mostly reside reside in the medium and 
deep sites. This indicates that during irradiation at 10K the HD molecules have sufficient 
mobility to hop from the shallow sites to deeper sites before they desorb. The similarity 
between the TPD curves obtained with HD and H+D irradiations indicates that molecules 
just formed on the surface quickly thermalize with the surface temperature and occupy the 
same adsorption sites as those molecules deposited from the gas phase [411 ] . 

In Fig. [3] we present several TPD traces obtained after irradiation with HD molecules at 
a sample temperature of 10K. The irradiation times are 30 (triangles), 60 (x), 120 (circles) 
and 240 (squares) s. In all cases the main peak is at 22K, as expected for irradiation with 
molecules, leading to first order kinetics. The high temperature peak appears around 31-32K 
for all exposures. 

In Fig. H] we present several TPD traces obtained after irradiation with H+D atoms at a 
sample temperature of 10K. The irradiation times are 30 (triangles), 60 (x), 120 (circles) and 
240 (squares) s. The location of the peak is clearly the same for all exposures, indicating a 
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first order kinetics. This peak temperature (around 23K) is only slighly higher than the 22K 
peak obtained for HD irradiation. However, the peak obtained for H+D irradiation exhibits 
a broader high temperature wing, which merges with the second peak. The first order 
feature of the H+D peak for irradiation at of 10K may indicate that some of the molecules 
are formed during irradiation or at early stages of the TPD run and quickly equilibrate 
with the surface. However, the broad high temperature wing may indicate that some other 
molecules are formed at later stages of the TPD run. In this case, the high temperature peak 
(which is difficult to resolve) may exhibit a combination of first and second order features. 
The experimental data was analyzed using the rate equation models described below. The 
results for the energy barriers for diffusion and desorption are summarized in Table I. 

IV. THE RATE EQUATION MODEL 

Consider an experiment in which beams of H and D atoms are irradiated on a surface. 
Atoms that stick to the surface hop as random walkers. The hopping atoms may either 
encounter each other and form H 2 , HD and D 2 molecules, or desorb from the surface in 
atomic form. As the sample temperature is raised, both the hopping and desorption rates 
quickly increase. The fluxes of H and D are denoted by /h and /d, respectively [in units of 
monolayer (ML) s -1 ]. These fluxes are related to the beam intensities by /h = F^/s and 
/d = Fv/s, where s (cm -2 ) is the density of adsorption sites on the surface. 

The process of molecular hydrogen formation is described by a rate equation model, 
which includes three surface reactions, namely H + H — > H 2 , H + D — > HD and D + D — > D 2 . 
The surface coverages of adsorbed H and D atoms are denoted by n H (ML) and riv (ML), 
respectively. Similarly, the coverages of the adsorbed H 2 , D 2 and HD molecules are denoted 
by riu 2 , ud 2 and n H D (ML). The time derivatives of the coverages of adsorbed atoms are 
given by the rate equations 



dnu 



fn - W H n H - 2a H nn - (a H + a D )n H n D 



dt 
driY) 



/d - WbUb - 2a D n5 - (a H + a jD )n H n D . 



(1) 



dt 



The time derivatives of the coverages of adsorbed molecules are given by 
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dn n 2 2 rrr 

—— = a H n H - Wu 2 n H2 

dn D2 2 

— — = a D n D - W D2 n D2 

dn UD 



(a R + a D )n H n D - W H d"hd- (2) 

The first terms on the right hand sides of Eqs. represent the incoming flux. For simplicity, 
we ignore the Langmuir rejection. The second terms in Eqs. (CQ) represent the desorption of 
H and D atoms from the surface, while the second terms in Eqs. (GJ) represent the desorption 
of molecules. The desorption coefficients are 

W x = vex P (-E*? s /k B T) (3) 

where v is the attempt rate (standardly taken to be 10 12 s~ l ), E x es is the energy barrier 
for desorption of species X, where X= H, D, H2, HD or D2, and T (K) is the surface 
temperature. The third and fourth terms in Eqs. ([I]) and the first terms in Eqs. (EJ) account 
for the formation of H2, HD and D2 molecules, where 

a x = vexp(-E x iS /k B T) (4) 

is the hopping rate of atoms of species X = H, D between adsorption sites on the surface 
and E x iS is the energy barrier for hopping. Here we assume that there is no energy barrier 
for the formation of hydrogen molecules on the surface, namely two hydrogen atoms that 
encounter each other form a molecule. The desorption rate Rx (ML s _1 ) of molecular species 
X is Rx = Wxnx- In particular, the rate of HD desorption is given by 

R HD = Whd^hd- (5) 

To simplify the analysis and reduce the number of fitting parameters it is desireable to 
use a model which includes only one isotope of hydrogen. Here we examine the effect of 
the isotopic difference on the TPD results. Due to its higher mass, the zero point energy 
for adsorbed D atoms is lower than for H atoms. As a result, the energy barriers for 
diffusion and desorption of D atoms are expected to be higher than for H atoms. The 
evaluation of the difference in the activation energies between the H and D isotopes requires 
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detailed knowledge of the atom-surface potentials. Such potentials are not available for the 
amorphous surfaces of interest here. However, one can obtain a rough idea about the isotope 
effect by considering the difference between the activation energies of H and D atoms on a 
graphite surface. In this case, measurements show that the laterally averaged binding energy 
on the basal plane of graphite, in the ground state level is 31.6 meV for H atoms and it is 
3.8 meV higher for D atoms 43j. The isotopic difference in the energy barriers for diffusion 
is expected to be at most as large as the difference in the energy barriers for desorption. 

Here we examine the effect of increasing the energy barriers for D diffusion and desorption, 
keeping those of H atoms unchanged. In Fig. [5] we present the TPD curves obtained from 
the rate equation model, following irradiation at 5.6K and heating at a constant rate of 
0.5K/s. Four curves are shown: the TPD results obtained when the parameters of H and 
D are identical, and given by Table I (solid line); the energy barrier for desorption of D 
atoms is raised by 4.0 meV (dashed-dotted line); the energy barriers for desorption and 
diffusion of D atoms are raised by 4.0 meV (circles); and the energy barriers for desorption 
and diffusion of both H and D are raised by 4.0 meV (dashed line). We find that raising the 
energy barriers for desorption (and diffusion) of D atoms gives rise to only a slight shift of 
the peak to a higher temperature. We conclude that the TPD curves are mainly determined 
by the isotope that interacts more weakly with the surface. It is thus appropriate to use a 
simplified model which includes only one isotope. The energy barriers obtained should be 
interpreted as the barriers for H atoms. 

In the models we assume a given density of adsorption sites on the surface. In terms of 
the adsorption of H atoms, all the adsorption sites are assumed to be identical, where the 
energy barrier for H diffusion is E^ s and the barrier for desorption is E^ s . This assumption 
is justified by the fact that such a simple model provides good fits to the TPD curves obtained 
after H+D irradiation. If there is a distribution of diffusion energy barriers, the hydrogen 
atoms are expected to be trapped most of the time in the deep adsorption sites. Thus, the 
diffusion-induced recombination is expected to be dominated by the deepest wells. This 
means that the parameters obtained from fitting the TPD curves to the rate equation model 
characterize the upper edge of the distribution of energy barriers. A more complete model 
that provides a connection between the roughness and the distribution of energy barriers for 



H diffusion and desorption was studied in Ref. 



44j. 



As for the adsorption of hydrogen molecules, we assume that the adsorption sites may 
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differ from each other. The energy barriers for desorption of HD molecules can be obtained 
directly from the TPD curves that follow HD irradiation. These curves, obtained for different 
surface temperatures during irradiation, reveal a total of three peaks, indicating that there 
are three types of adsorption sites. In general, we assume that the adsorption sites are 
divided into J types according to the binding energies of trapped HD molecules. Ignoring 
the isotopic differences for molecules, the energy barrier for desorption of H2 molecules from 
an adsorption site of type j is E#®(j), where j — 1, . . . , J. The parameter fij represents the 
fraction of the molecules that are trapped in sites of type j upon formation or adsorption and 
Y^j (J>j = 1- Let nn 2 0) (ML) be the coverage of H 2 molecules that are trapped in adsorption 
sites of type j, where j = 1, . . . , J. The rate equation model takes the form [39j 



nn = /h (1 - n u ) - Wnn u - 2a H n H 2 (6a) 
nn 2 (j) = Hja u n H 2 -W U2 (j)n H2 (j), (66) 

where j = 1, . . . , J. The first term on the right hand side of Eq. ( !6al) represents the incoming 
flux in the Langmuir kinetics. In this scheme H atoms deposited on top of H atoms already 
on the surface are rejected. There are indications that adsorbed H 2 molecules do not lead 
to Langmuir rejection of H atoms [3]. The parameter /h represents an effective flux (in 
units of ML s _1 ), namely it already includes the possibility of a temperature dependent 
sticking coefficient. The second term in Eq. fl6a[) represents the desorption of H atoms from 
the surface. The third term in Eq. (16aj) accounts for the depletion of the H atoms on the 
surface due to diffusion-mediated recombination into H 2 molecules. Eq. (16fo]) accounts for 
the population of molecules on the surface. The first term on the right hand side represents 
the formation of H 2 molecules that become adsorbed in a site of type j. The second term 
in Eq. ([65]) describes the desorption of H 2 molecules from sites of type j, where 

W H2 (j) = vexp[-E^(j)/k B T} (7) 
is the H 2 desorption coefficient. The H 2 production rate _Rh 2 (ML s _1 ) is given by: 

j 

^ Ha = X) W H 2 0>H 2 0')- ( 8 ) 
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To analyze the diffusion of H2 molecules, one needs to use a more complete model 39( . This 
model includes additional parameters, namely the energy barriers E^(j) for hopping of H2 
molecules out of sites of type j, as well as the partial densities Sj (cm -2 ) of such sites. It is 
difficult to extract, from the experimental data, unique values for all these parameters with 
sufficient confidence. 



V. ANALYSIS OF THE EXPERIMENTAL RESULTS 



In Ref. (4l| we presented a series of TPD traces obtained after irradiation of H and 
D atoms on an amorphous silicate sample at low surface temperatures (5.6K). Each trace 
exhibits a large peak at a lower temperature and a small peak at a higher temperature. The 
location of the low temperature peak shifts to the right as the irradiation time is reduced, 
suggesting second order kinetics. We should note that a somewhat similar behavior was 
bund in TPD experiments in which D2 molecules were irradiated on amorphous water ice 



321 . l33| . In that case the second-order like behavior was interpreted as a result of the 
saturation of the deepest adsorption sites. However, in our experiments the coverage of 
hydrogen atoms on the surface is very low. Therefore, it is unlikely that saturation effects 
play a significant role. 

This observation indicates that the molecules are formed during the heating and not 
during the irradiation stage. Thus, the adsorbed hydrogen atoms are immobile during ir- 
radiation at 5.6K. This means that tunneling alone is not sufficient in order to provide 
significant mobility to the adsorbed hydrogen atoms, and their mobility is dominated by 
thermal activation. A similar conclusion was reached from analysis of data from polycrys- 



talline olivine 



36). 



Assuming that the Langmuir rejection mechanism applies, hydrogen atoms that are de- 
posited on top of already adsorbed atoms are rejected. This provides a prediction for the 
coverage of adsorbed atoms after irradiation time t^. Taking the Langmuir rejection into 



account, the coverage is given by [45] 



^h(*o) = 1 - exp(-/ H • to)- (9) 

In the experiment, the coverage after irradiation can be evaluated using the total yield of 
HD molecules in each TPD run. The total yields for irradiation at 5.6K and several exposure 
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times were evaluated and fitted according to Eq. (Q. It was found that the flux of incoming 
atoms is f H = 7.0 x 10" 4 (in ML s -1 ). Since the beam intensities Fn and i*D are known, 
one can use the relations fx = F x /s, where X=H, D, to obtain the density of adsorption 
sites. It was found that the density of adsorption sites on the amorphous silicate sample is 
s = 7 x 10 14 (sites cm -2 ). 

The experimental results were fitted using the rate equation model described above. 
The parameters for the diffusion and desorption of hydrogen atoms and molecules on the 
amorphous silicate surface were obtained. These include the energy barrier E^ = 35 (meV) 
for the diffusion of H atoms and the barrier E^ s = 44 (meV) for their desorption. The value 
obtained for the energy barrier for desorption should be considered only as a lower bound, 
because the TPD results are insensitive to variations in E^ s , as long as it is higher than the 
reported value. The desorption energy barriers of HD molecules adsorbed in shallow (lower 
temperature peak) and deep (higher temperature peak) sites, are given by E^(l) = 35 and 
E^(2) = 53 (meV). 

The rate equation model is integrated using a Runge Kutta stepper. For any given choice 
of the parameters, one obtains a set of TPD curves for the different irradiation times used in 
the experiments. The actual time dependence of the temperature, T(t), recorded during the 
experiment, is taken into account as follows. The rate equations [Eqs. (6)] are integrated 
numerically using a Runge Kutta stepper. The rate constants for diffusion and desorption 
[given by Eqs. (3), (4) and (7)], which depend on the temperature, are adjusted during 
the integration according to the experimentally recorded temperature curve, T(t). Thus the 
simulation fully reflects the physical conditions during the experiment. 

In the first step, the barriers E^(j), j = 1, . . . , J, for the desorption of molecules are 
obtained using the results of the experiments in which HD molecules are irradiated on the 
surface. To obtain better fits of the peak shape, suitable Gaussian distributions of e nerg y 
barriers around these two values can be incorporated (see Figs. 1 and 2 in Ref. 41|). 
However, even with these Gaussian distributions it is difficult to fit the leading edge of the 
TPD curves. This may indicate that the distribution of the energy barriers is non-symmetric 
and exhibits a broader tail on the lower side. 

A wide range of energies can be expected for a morphologically heterogeneous surface. For 
example, in a calculation of H adsorption energies and diffusion barriers on a water cluster 
simulating the surface of amorphous water ice, Buch and co-workers obtained Gaussian 
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distribution for these energies 




Using a single value of each barrier, one obtains fits which consist of rather narrow peaks. 
Their locations coincide with the experimental peaks but they do not capture the tails. 

The results of the H+D experiments (circles in Fig. 1) are fitted using the rate equation 
model (solid line) with the activation energies specified above. The weights of the molecular 
adsorption sites are chosen to be \L\ = 0.6 and ^2 = 0.4. These values reflect the total 
relative yields of desorbed HD molecules in the two peaks. The results of H+D irradiation 
at 10K, shown in Fig. [2j were also fitted using the rate equation model. The activation 
energies, E# and E^ s are the same as above. However, in this case it is found that the 
molecular adsorption sites indexed by j = 1 are too shallow to trap molecular hydrogen when 
the irradiation is done at 10K. To account for this fact we impose \x\ = 0. The TPD curve 
is fitted with E$*(2) = 53 (as above) and E^(3) = 75 (meV). The weights are ^2 = 0.5 
and /z 3 = 0.5. 



VI. DISCUSSION AND APPLICATIONS 



The processes of molecular hydrogen formation on amorphous silicate surfaces at low 
temperatures are based on physisorption forces. The activation energies that we obtained 
show that these values are higher than expected for other physisor ptio n interactions that 
have been measured between H and the surface of single crystals [48j. This is probably 
due to the more complex chemical and morphological composition of the amorphous silicate 
samples. Using the parameters obtained from the experiments we now calculate the recom- 
bination efficiency of hydrogen on amorphous silicate surfaces under interstellar conditions. 
The recombination efficiency is defined as the fraction of the adsorbed hydrogen atoms which 
come out as molecules, namely 49| 



n 



R 



H 2 



C/h/2)- 

Under steady state conditions, the recombination rate Ru 2 (ML s _1 ) is given by i? H 
oh^h 2 , where 



(10) 



4a H 



-(W H + / H ) + V(W U + / H ) 2 + 8a H / H 



(11) 
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Within the assumptions of this model, the binding energies E#* of H2 molecules to the 
surface do not affect the rate of H2 formation. They only determine the residence time of 
H 2 molecules on the grain and thus affect the steady state coverage of H 2 molecules on the 
surface. 

In Figs. El^a) and MJo) we present the recombination efficiency and the coverage of H 
atoms, respectively, vs. surface temperature for an amorphous silicate sample under flux of 
/ H = 5.2 x 1(T 9 (ML s" 1 ). This flux is within the typical range for diffuse interstellar clouds, 
where bare amorphous silicate grains are expected to play a crucial role in H 2 formation. 
This flux corresponds to gas density of 100 (atoms cm -3 ), gas temperature of 100K and 
a density of 7 x 10 14 adsorption sites per cm 2 on the grain surface. A window of high 
recombination efficiency is found between 9-14K, compared to 6-10K for polycrystalline 
silicate under similar conditions. The surface coverage of H atoms decreases dramatically 
within the efficiency window, from nearly full monolayer at 9K to about 10~ 4 ML at 14K. 
This is due to the fact that as the temperature increases, the surface mobility of H atoms is 
enhanced and their residence time before recombination is reduced. 

In g eneral, the high efficiency window for hydrogen recombination is bounded from below 

by y 



E diS 
k B (\nu - ln/ H ) 



and from above by 



2^des ^diff 

^hi g h(/H) = "TTi i — FT' ( 13 ) 

k B (\nv- ln/ H ) 

These bounds are obtained by solving Eq. (l6a|) under steady state conditions for recombina- 
tion efficiency of 50%. Thus, inside the window the efficiency is higher than 50%, while out- 
side it is lower than 50%. The recombination efficiency declines sharply near these bounds. 
The width of the high efficiency window is proportional to the difference, E^ s — E^ 5 , be- 
tween the diffusion and desorption energy barriers of H atoms. The location of this window 
exhibits a logarithmic dependence on the flux /h and slowly shifts to higher temperatures 
as /h increases. 

At temperatures higher than T^gh atoms desorb from the surface before they have suffi- 
cient time to encounter each other. At temperatures lower than Ti ow diffusion is suppressed 



16 



and the surface is saturated by immobile H atoms. As a result, the Langmuir-Hinshelwood 
mechanism is no longer efficient. Saturation of the surface with immobile H atoms might 



render the Eley-Rideal mechanism more efficient in producing some recombination [36l. 139]. 
Our results thus indicate that the recombination efficiency of hydrogen on amorphous sil- 
icates is high in a temperature range, which is relevant to interstellar clouds. Therefore, 
amorphous silicates seem to be good candidates for interstellar grain components on which 
hydrogen recombines with high efficiency. 

To apply the results of the experiments to H 2 formation in the ISM, consider for simplicity 
a spherical grain of radius r. The cross section of such grain is o = nr 2 and the number 
of adsorption sites on its surface is S — 4irr 2 s. We denote by N-r the number of H atoms 
adsorbed on such grain. Its time dependence is given by 

N n = F R a - W R N n - A n N^ (14) 

The incoming flux density is given by Fh = n^ as t>H, where n^ as (cm -3 ) is the density of H 
atoms in the gas phase and vn is their average velocity. The parameter A-r = an/S, is the 
sweeping rate, which is approximately the inverse of the time it takes for an H atom to 
visit nearly all the adsorption sites on the grain surface (for a more precise evaluation of the 



sweeping rate see Refs. 



51 



52|). 



The production rate of H 2 molecules on a single grain, is given by 

R& 2 = A H N R 2 . (15) 

Under given flux and surface temperature, for grains that are large enough to hold many 
H atoms, Eqs. provide a good description of the recombination process. However, in 
the limit of small grains and low flux, Nn may be reduced to order unity or less. Under 
these conditions, Eqs. ([6]) become unsuitable, because they ignore the discrete nature of the 
population of adsorbed atoms and its fluctuations. 

To account for the reaction rates on small grains, a modified set of rate equations was 



introduced and applied to a variety of chemical reactions [53j, |54j. In these equations the 
rate coefficients are modified in a semi-empirical way taking into account the effect of the 
finite grain size on the recombination process. The modified rate equations take into ac- 



count correctly the length scales involved in the recombination process on small grains 



50] 
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However, they still involve only the average values and ignore fluctuations, thus providing 
only approximate results for the reaction rates on small grains. 

To account correctly for the effects of fluctuations on the recombination rate on small 
grains, simulations using the master equation are required. The dynamical variables of the 
master equation are the probabilities P(Ng) of having a population of Ah hydrogen atoms 
on the grain. In the case of hydrogen recombination, the master equation takes the form 



45|, 



P(N K ) = F H [P(N U - 1) - P(N U )] + W u [(N H + 1)P(N U + 1) - N R P(N R )] 

+ A H [(N R + 2)(iV H + l)P(N u + 2) - N R (N U - l)P(N K )], (16) 

where Ah = 0, 1, 2, ... , A^f ax (the master equation must be truncated in order to keep the 
number of equations finite). The first term on the right hand side of Eq. fflBT) describes the 
effect of the incoming flux. The probability P(A^h) increases when an H atom is adsorbed by 
a grain that already has Ah — 1 adsorbed H atoms, and decreases when it is adsorbed on a 
grain with A^h atoms. The second term accounts for the desorption process. The third term 
describes the recombination process. The recombination rate is proportional to the number 
of pairs of H atoms on the grain, namely Ah(Ah — l)/2. Therefore, the H2 production rate 
per grain can be expressed in terms of the first two moments of P(Nn), according to 

K 2 =M(N 2 H )-(N R )), (17) 

where 

«> = E N H k P(N R ). (18) 

N H =0 

The master equation can be simulated either by direct numerical integration or by a 



stochastic implementation via Monte Carlo (MC) methods [56l . |57| | . A significant advantage 
of the direct integration over the MC approach is that the equations can be easily coupled 
to the rate equations of gas-phase chemistery. However, the number of coupled equations 
increases exponentially with the number of reactive species, making direct integration infea- 



sible for complex reaction networks of multiple species 



58, 



591 ] . To solve this problem, two 
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approaches have been proposed: the multiplane method 60|, |6l| and the moment equations 



62. 



63 



64]. 



The multiplane method 60| is based on breaking the network into a set of maximal 
fully connected subnetworks (maximal cliques). It involves an approximation, in which 
the correlations between pairs of species that react with each other are maintained, while 
the correlations between non-reacting pairs are neglected 61j. The result is a set of lower 
dimensional master equations, one for each clique, with suitable couplings between them. 
For sparse networks, the cliques are typically small and mostly consist of two or three nodes. 
This method thus enables the simulation of large networks much beyond the point where 
the master equation becomes infeasible. 

The moment equations include one equation for the population size of each reactive 
species (represented by a first moment) and one equation for each reaction rate (represented 
by a second moment). These equations are obtained by taking the time derivative of each 
moment and using the master equation to express the time derivatives of the probabilities 
62|. In the resulting equations, the time derivative of each moment can be expressed as a 
linear combination of first, second and third order moments. To close the set of moment 
equations one must express the third order moments in terms of first and second order 
moments. This is achieved by the incorporatin of a suitable truncation scheme of the master 



equation 



63 



641 ] . The moment equations provide the most efficient incorporation of 



stochastic grain chemistry, including hydrogen recombination and other surface reactions 
into models of interstellar chemistry. 



VII. SUMMARY 



We have reported the results of experiments on molecular hydrogen formation on amor- 
phous silicate surfaces, using temperature-programmed desorption. In these experiments 
beams of H and D atoms were irradiated on the surface of an amorphous silicate sample. 
The desorption rate of HD molecules was monitored using a mass spectrometer during a 
subsequent TPD run. The results were analyzed using rate equations. It was shown that a 
model based on a single isotope provides the correct results for the activation energies for 
diffusion and desorption of H atoms. The barriers E^ s , E^ s and E^, as well as the density 
of adsorption sites on the surface were obtained. 
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The results were used in order to evaluate the formation rate of H 2 on dust grains under 
the actual conditions present in interstellar clouds. It was found that under typical conditions 
in diffuse interstellar clouds, amorphous silicate grains are efficient catalysts of H 2 formation 
when the grain temperatures are between 9-I4K. This temperature window is within the 
typical range of grain temperatures in diffuse clouds. It was thus concluded that amorphous 
silicates are efficient catalysts of H 2 formation in diffuse clouds. The recently developed 
computational methodologies for the evaluation of reaction rates on interstellar grains were 
briefly reviewed. 
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TABLE I: Parameters for molecular hydrogen formation: the diffusion and desorption barriers of 
H atoms, the desorption barriers of HD molecules and the density of adsorption sites. 

Material ^ iff (meV) ^ cs (meV) Egg(meV) s (sites cm" 2 ) 

Amorphous Silicate 35 44 35,53,75 7 xlO 14 
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FIG. 1: TPD curves of HD desorption after irradiation H+D atoms (circles) and HD molecules 
(x) on the amorphous silicate sample at surface temperature of 5.6K. The fit of the TPD curve for 
H+D irradiation, obtained using the rate equations, is also shown (solid line). The temperature 
ramp T{t) during the TPD stage is shown in the inset. 
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FIG. 2: TPD curves of HD desorption after irradiation H+D atoms (circles) and HD molecules 
(x) on the amorphous silicate sample at surface temperature of 10K. The fit of the TPD curve for 
H+D irradiation, obtained using the rate equations, is also shown (solid line). The temperature 
ramp T{t) during the TPD stage is shown in the inset. 
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FIG. 3: TPD traces obtained after irradiation with HD molecules at a sample temperature 
10K. The irradiation times are 30 (triangles), 60 (x), 120 (circles) and 240 (squares) seconds. 
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FIG. 4: TPD traces obtained after irradiation with H+D atoms at a sample temperature of 10K. 
The irradiation times are 30 (triangles), 60 (x), 120 (circles) and 240 (squares) seconds. 
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FIG. 5: Simulated TPD curves of HD production and desorption obtained from Eqs. ([I]), following 
irradiation by H+D atoms at 5.6K and heating at a constant rate of 0.5K/s. In this Figure it is 
assumed that HD molecules desorb upon formation. Four choices of parameters are shown: the 
parameters of H and D are identical, and are given by Table I (solid line); the desorption energy 
barrier of D is raised by 4.0 meV (dashed-dotted line); both the desorption and diffusion energy 
barriers of D are raised by 4.0 meV (circles); the desorption and diffusion energy barriers of both 
H and D are raised by 4.0 meV (dashed line). The results indicate that the higher binding energy 
of D atoms has little effect on the TPD curves. The peak heights are normalized to 1, to make it 
easier to compare the peak locations. 
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FIG. 6: (a) Calculated recombination efficiency and surface coverage (b) of hydrogen on the 
amorphous silicate sample vs. surface temperature, under steady-state conditions. A temperature 
window of high efficiency is found. Within this window, the surface coverage is reduced from a full 
monolayer to about 10~ 4 ML. 
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